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It is shown that a local bend of a nanowire is a source of pinning potential for a transversal 
head-to-head (tail-to-tail) domain wall. Eigenfrequency of the domain wall free oscillations at the 
pinning potential and the effective friction are determined as functions of the curvature and domain 
wall width. The pinning potential originates from the effective curvature induced Dzyaloshinsky- 
like term in the exchange energy. The theoretical results are verified by means of micromagnetic 
simulations for the case of parabolic shape of the wire bend. 

PACS numbers: 75.10.Hk, 75.10.Pq, 75.40.Mg, 75.60.Ch, 75.78.Cd, 75.78.Fg 


I. INTRODUCTION 

Dynamics of domain walls in nanowires underlies mod¬ 
ern memory^’^ and logic^^® devices, see also review Ref. 6. 
Most of these designs require wires with curvilinear seg¬ 
ments, e.g. vertical configuration of the racetrack mag¬ 
netic memory.^ Therefore, the modification of dynamic 
and static properties of a domain wall due to the wire 
curvature is of high importance for the applications. It is 
known that the local curvilinear defects of the wire can 
be a source of pinning potential for the domain wall.^^® 
However, in contrast to widely studied mechanisms of the 
domain wall pinning, such as interaction with artificial 
notches,surface roughnesses,inhomogeneities 
in the magnetocrystalline anisotropy distribution, 
the role of the curvature remains unclear and the gen¬ 
eral theory is absent. Nevertheless, one should note some 
progress in studying the simplest case of circular wire ge¬ 
ometry, where the domain wall pinning is achieved due 
to the applied external magnetic field. 

As it was recently shown^® the exchange interac¬ 
tion in a curvilinear wire can produce an effective 
Dzyaloshinsky-like term. In this paper we demonstrate 
that the inhomogeneous distribution of this curvature in¬ 
duced Dzyaloshinsky-like contribution can be an origin of 
domain wall pinning at the local wire bend, analogously 
to the case of studied anisotropy inhomogeneities. 
Considering the coordinate dependent curvature of the 
wire we propose a general approach valid for a wide class 
of geometries. 


II. MODEL AND GENERAL RESULTS 

In this paper we consider a plane curved ferromagnetic 
wire of circular cross-section. Such a wire can be param¬ 
eterized in the following way: 

r{s,x,p) = Tis) + pcosxe^is) + psinxeBis). (1) 

Here the three-dimensional radius vector r defines the 
space domain, occupied by the wire, the two-dimensional 


vector 7(s) = 'yxis)x+Xy{s)y determines the wire central 
line, which lies within the xy-plane, with s being the nat¬ 
ural parameter (arc length). Within a wire cross-section 
the polar coordinates p S [0, R] and x € [0, 27r) are used, 
where R is the wire radius. We use here the Frenet-Serret 
basis (bt, 6^,63) with 63 = 7'(s), = -y”{s)/ k{s), and 

63 = 63 X gpj being the tangential, normal, and binor¬ 
mal unit vectors, respectively. Here the wire curvature 
k{s) = |7"(s)| is introduced. 

Using the Frenet-Serret basis one can introduce the 
angular magnetization parameterization 


m = cos 6 e-j. + sin 0 cos ^6^-1- sin 0 sin </> 63, (2) 


where m = M/Mg is the normalized magnetization unit 
vector with Mg being the saturation magnetization. 

Magnetization dynamics of this system can be studied 
by means of phenomenological Landau-Lifshitz equations 


— sin 06 = ujq— a sin^ i 

d(j) 


(5£ 

sm0(j) = (jjo^ + a0, (3) 
do 


where overdot indicates the time derivative, frequency 
Wo = 47r7oMs determines the characteristic time scale 
of the system with 70 being the gyromagnetic ratio, 
£ = Ej^nMl is normalized total energy, and a is the 
damping coefficient. 

We start with a simple model, which takes into account 
only two contributions to the total magnetic energy 


+00 

£=§ 

— 00 


fct cos^ 0] ds, 


( 4 ) 


namely, exchange one and easy-tangential anisotropy 
- the second term in (4). Here § = Tri?^ is area of the wire 
cross-section, £ = ^A is the exchange length with 
A being exchange constant, and fct = K / -1-1/4 
is the dimensionless anisotropy constant. Here K > 
0 is constant magnetocrystalline anisotropy of easy- 
tangential type and the term 1/4 comes from the mag¬ 
netostatic contribution. It is known^^^^° that for thin 
wires of circular and square cross-sections the magneto¬ 
static energy is reduced to an effective easy-tangential 
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FIG. 1. (Color online) Equilibrium state of the transversal do¬ 
main wall at the wire bend. Top row shows the magnetization 
distribution obtained using micromagnetic NMAG simulations 
for permalloy parabolic wire with R = 5 nm, total length 
L — 1 fim and kq = 0.05 nm“^ (only the bend vicinity is 
shown). The bottom row demonstrates comparison of mag¬ 
netization components mt = m e^ and rUn = m et, obtained 
from the simulations (markers) and from the Ansatz (7) 
(lines). In the latter case the domain wall width A = 2£ 
was used. 


shape anisotropy with Kf.s = including case of 

a curvilinear wire.^® Competition of the exchange and 
anisotropy contributions results in the length scale of the 
system w = For magnetically soft wires K = 0 

and w = 2£. In (4) and everywhere below we restrict 
ourselves with the case of thin wire R w, thus we as¬ 
sume that the magnetization varies along the wire and it 
is uniform within a wire cross-section: m = m{s). 

In terms of the angular parametrization (2) the ex¬ 
change energy density has a form^® 

<^ex = (O' + K cos (/))^ -I- {(j)' sin 6* — K cos 9 sin c/))^ , (5) 

where the prime denotes derivative with respect to s. In 
(5) it is taken into account that a plane wire has zero 
torsion. 

Let us first analyze static solutions of (3). In this case 
there is a solution (j) = (j)^ = 0, ir, which corresponds 
to a planar magnetization distribution: vectors m(s) lie 
within the wire plane. The corresponding function 9 is 
determined by an inhomogeneous Sine-Gordon equation 
(for details see Appendix A) 

w^O" — sinOcosO =—w^n'cos(j)o. (6) 

For the case k' = 0 (rectilinear or circular wire) 
Eq. (6) has the well known domain wall solutions cos 9 — 
—ptanh[(s — q)/A] of head-to-head (p = 1) or tail-to-tail 
(p = —I) types. Here q determines the domain wall po¬ 
sition and A = w is the wall width. For the case k' ^ 0 
an additional driving force appears and the curvature 


induced domain wall dynamics is expected. In the fol¬ 
lowing, we consider a case of a localized curvature, i.e. 
«:(±oo) = 0 and k'(±oo) = 0, see Figs. I, 2(a), this re¬ 
sults in the boundary conditions cos6>(±oo) = ^p, the 
same as for a rectilinear wire. We also restrict ourselves 
with the case w'^k' I assuming that the domain wall 
remains its form and the curvature influence results in a 
weak modihcation of the width A. Therefore, to analyze 
the domain wall properties (static as well as dynamic 
ones) we use collective variable approach®^’®^ based on 
the simple g — $ model®^’®® 


cos 6* = 


—p tanh 


s - qjt) 
A 


(/) = $(t). 


( 7 ) 


which is widely used for different types of domain walls 
and various drivings.The domain wall po¬ 
sition q and phase $, which determines orientation of 
the transversal magnetization component, make a canon¬ 
ically conjugated pair of collective variables. The domain 
wall width A is shown to be a slaved variable,i.e. 
A{t) — A[g(t), <i)(t)]. However, this is not the case when 
$ depends on coordinated^ (it is possible for a three di¬ 
mensional wire, when the torsion is present), in this case 
a generalization of the g — $ modeled must be applied. 

Substituting now the Ansatz (7) into (4) one obtains 
the total energy in the form (up to an additive constant) 

-I- ni'^p cos $k(< 7 ) — ^®Asin® $«:^(g), 

( 8 ) 

where the condition kA <C I was applied when in¬ 
tegrating (4), for details see Appendix B. The first 
term in (8) reflects the exchange-anisotropy competition 
which determines the domain wall width for a rectilin¬ 
ear wire. The second and the third terms originate 
from Dzyaloshinsky-like and anisotropy-like terms, re¬ 
spectively, which effectively appear in the exchange en¬ 
ergy due to the curvature.®® Minimization of (8) with 
respect to parameters q and $ results in the following 
equilibrium values qg and 

K'iqo)=0, cos$o = -p. (9) 



Thus, the domain wall is pinned at the maximum of 
the curvature, and the phase selectivity takes place: 
the head-to-head domain wall is always directed out¬ 
ward while the tail-to-tail one is directed inward the 
bend, see Fig. 1. The latter effect was recently observed 
experimentally.^^ There is an intuitive explanation: the 
choice of $0 defined by (9) makes the magnetization dis¬ 
tribution more homogeneous, see Fig. 1, which decreases 
the exchange energy. 

As it follows from (8), the equilibrium value of the 
domain wall width Aq — w is the same as for a rec¬ 
tilinear wire. However, if the wall is perturbed and 
the values of quantities (q, $) deviate from the equilib¬ 
rium (9), the domain wall width is modihed as follows 
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A{q,^) = £[kt — fcb(( 7 ) sin^ where fcb = pK^{q) 

is coefficient of the curvature induced effective easy- 
binormal anisotropy. The similar domain wall width 
modification caused by a geometrical constrain was dis¬ 
cussed in Ref. 42. 

As it follows from the Fig. 1 the used Ansatz (7) and 
the obtained static results (9) are in a good agreement 
with the simulations results. 

Let us now proceed to dynamical properties of the do¬ 
main wall. In terms of the collective variables the equa¬ 
tions of motion (3) take a form (see Appendix B) 


Wo 9£ . • 


■ Wo 9£ a . , 


Our goal is to study linear dynamics of the domain wall 
in vicinity of the equilibrium position. With this purpose 
we introduce small deviations in the way q{t) = qo + q{t) 
and = 4)0 -I- l>(t). For the limit case kAq —)• 0 the 
equations of motion ( 10 ) linearized with respect to the 
deviations read 


In this case the equilibrium position q^ is defined by the 
equation 


+ 00 


«'(s) 


cosh 


ds = 0 


(15) 


and the equilibrium value of $0 coincides with one de¬ 
fined in (9), for details see Appendix B. 

It is important to emphasize that the pinning poten¬ 
tial and the corresponding domain wall oscillations have 
purely exchange origin, while the weak interactions, e.g. 
magnetostatic one, contribute mediately via value of pa¬ 
rameter Ag. 


III. THE CASE OF PARABOLIC BEND 

Let us now check the obtained general results for an 
example of parabolic wire, whose central line has a form 


(l + a^) 


q 

$ 


Wg^^TT 


0 

k"(9o) 




—a 


Aq 



( 11 ) 


For the case of small a the solution of (11) results in 
harmonic decaying oscillations q oc sin(Ot-|-i5g)e“''*, l> oc 
cos(r2t -I- 6o)e~'^^ with frequency 


n « Wg£^7ri/)t(gg)|K"((7g)| (12) 


and modified effective friction 




i‘^K{qo) 


7T t K 


(13) 


The phase ^g is determined by the initial conditions. In 
(12) and (13) it is taken into account that K(gg) > 0 and 
k"{ qo) < 0 at the maximum point gg. 

It should be noted that an energy expression analogous 
to ( 8 ) was earlier obtained for a certain case of a circu¬ 
lar wire segment.However, the circular geometry does 
not produce any geometrical pinning potential due to the 
constant curvature k = const, which results in H = 0 and 
undefined gg. Nevertheless, the domain wall can have an 
equilibrium position in circular segment in the presence 
of an external magnetic field. 

In more general case AgK G (0,1) the Ansatz (7) leads 
to the following expressions for the eigen-frequency and 
effective friction 




A = A, 


A2 

+ 00 

2 f k"(s) 


A2 

^0 


— 00 
+00 


cosh 


-ds, 


® = 


k(s) 


+ 00 


cosh 


-ds — An 


K^{S) 


cosh^ 


ds. 


(14) 


y = XX 


"oyy- 


(16) 


The wire curvature k = «;g(l -I- k^x^) as a func¬ 
tion of natural parameter s = Kq^/{kqx), with f{^) = 

5 C-\/l + + arcsinh(^) , is shown in Fig. 2(a). It has 

a single maximum at point sg = 0 with extreme value 

K(sg) = Kg. 

In accordance with (15) the equilibrium domain wall 
position is gg = 0 , i.e. the domain wall is located at the 
extreme point of the bend. The equilibrium phase value 
4)g is determined by (9), see Fig. 1. 

In the narrow domain wall limit AgKg —> 0 the eigen- 
frequency ( 12 ) and effective friction read 


H « a;g-\/37r(Kg£)^, 77 « aa;g7r£Kg/4. (17) 


In order to verify our analytical results we performed 
two types of numerical simulation: (i) micromagnetic 
simulations of magnetically soft nanowire (NMAG code); 
(ii) numerical solution of the Landau-Lifshitz equations 
for a discrete chain of magnetic moments. 

Let us start by considering the magnetically soft wires, 
whose geometries are determined by (1) and (16). The 
radius of cross-section R = 5 nm and length L = 1 /rm are 
fixed for all studied samples, while the extreme curvature 
Kg is varied in the range Kg G [0.005, 0.05] nm“^ with 
step AKg = 0.005 nm“^. 

Permalloy is chosen as a material with the following 
parameters: exchange constant A = 13 pJ/m, satura¬ 
tion magnetization Mg = 860 kA/m, and damping coeffi¬ 
cient a = 0.01. These parameters results in the exchange 
length i « 3.7 nm and Wg = 30.3 GHz. Thermal effects 
and anisotropy are neglected. 

The magnetization dynamics is studied by means of 
numerical simulation of the Landau-Lifshitz equation ap¬ 
plying the NMAG code^^, wherein an irregular tetrahedral 
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FIG. 2. (Color online) (a) Curvature of the parabolic wire (16) with kq — 0.05 nm“^. Insets (b) and (c) shows the time 
dependences (with the corresponding Fourier spectra) of position g(t) and phase <E>(f) of the domain wall, perturbed at vicinity 
of the equilibrium point go = 0. The data are obtained by means of micromagnetic simulations for a wire with kq = 0.05 nm“^ 
and R = 5 nm. 


mesh with cell size about 1.75 nm is used. Only three 
magnetic interactions were taken into account, namely 
exchange, magnetostatic and Zeeman contributions. 

The numerical experiment consists of three steps. Ini¬ 
tially, we relax the domain wall structure in an over¬ 
damped regime (a = 0.1) in order to determine the 
equilibrium values of collective variables: position and 
phase 4)0. The obtained results fully coincides with the 
prediction (9). To determine the values of q and 4) we 
extract the curvilinear magnetization components nit = 
m-eT, TOn = m-eN, and Wb = m e^ from the simulation 
data, and apply fitting with the Ansatz (7). Namely, the 
position q is determined as a fitting parameter for the 
function mt(s) = —ptanh[(s — g)/A], then the phase is 
determined from the equation tan4) = mh{q)/rnn{q)■ 

In the second step we slightly perturb the domain wall 
phase 4> by applying a weak magnetic field B = BqEz 
perpendicularly to the wire plane, where Bq = 25 mT. 
After the system relaxation in the applied field B the 
field is switched off and the magnetization dynamics is 
simulated for the natural value of the damping coefficient 
(the third step). Since the variables 4) and q are canoni¬ 
cally conjugated, see Eqs. (10), the perturbed dynamics 
of 4) induces the dynamics of q, i.e. the domain wall 
starts to move. Typical time dependences of these col¬ 
lective variables are shown in Fig. 2(b,c), where one can 
see harmonic decaying oscillations with well pronounced 
frequency Q and friction parameter rj. In this way we 
determine and q for all range of the studied curva¬ 
tures Kq. The frequency If as well as the friction q in¬ 
creases with the curvature increasing, see Fig. 3 (round 
markers). One should note a good agreement of the ob¬ 
tained numerical results with the analytical prediction 
(14) with Aq = 2£, see solid lines in Fig. 3 and Ap¬ 
pendix B for details. This means that the approxima¬ 
tion of magnetostatic interaction by the effective easy- 
tangential anisotropy is physically sound for a domain 
wall dynamics in the thin nanowires. 



FIG. 3. (Color online) Eigenfrequency of the domain wall os¬ 
cillations in vicinity of the equilibrium - the extreme point of 
parabolic bend, see Fig. 1. Solid and dashed lines correspond 
to the predictions (14) and its narrow-wall asymptotics (17), 
respectively. Markers shows the results of numerical simula¬ 
tions for nanowires (disks) and discrete chains of magnetic 
moments (triangles). The inset demonstrates dependence of 
the effective friction on the curvature parameter. 

Additionally, we study dynamics of discrete chains of 
magnetic moments mi, with i = 1,A, aligned along 
parabolic line (16). The moments are aligned equidis- 
tantly with the fixed step size As along the chain. As pre¬ 
viously, only three magnetic interactions are taken into 
account, namely exchange (with £ = 3As), dipole-dipole 
and Zeeman contribution. Total number of moments is 
fixed N = 76 while the curvature parameter kq is varied 
but the restriction kqAs ^ 1. Dynamics of this system is 
described by a set of N vector Landau-Lifshitz equations 
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which a solved by means of common numerical methods, 
in details this procedure is described in Ref. 44. 

The domain wall oscillations are studied in the same 
way as for the wires described above. The resulting eigen- 
frequencies and friction parameters are shown in Fig. 3 
by triangles. 

In conclusion, we demonstrate an effect of curvature 
induced pinning of a transversal head-to-head (tail-to- 
tail) domain wall at a local wire bend. The curvature 
induced pinning potential has purely exchange origin. 
We obtain an expression (14) for frequency and effec¬ 
tive friction of the free domain wall oscillations within 
the pinning potential, and we demonstrate its validity 
numerically for the range of parameters kA ^ 10“^. 
For the case of a narrow domain wall or weak curva¬ 
ture {kA ~ 10“^) the approximations (12) and (13) can 
be used. The good agreement of analytical predictions 
and results of full scale numerical simulations for mag¬ 
netically soft cases {K = 0) demonstrates that the ap¬ 
proximation of magnetostatic interaction by the effective 
easy-tangential anisotropy is physically sound for a do¬ 
main wall dynamics in the thin nanowires. 
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Appendix A: Static solutions for a planar wire 


for Lagrange function^® 

OO 

C = —— f ^sin6*0ds—£ (B2) 

Wo j 

— OO 


and dissipative function 


OO 



— OO 


Substituting the Ansatz (7) into (B2) and (B3) and per¬ 
forming the integration over s (along the curve) one ob¬ 
tains the effective Lagrange and effective dissipative func¬ 
tions in the form 


r e rr _ ^ 

fC eff — o, J eff — CH 

Wo Wo 


rz2 


A 


Ad) 


k2 


(B4) 


Substituting (B4) into the Euler-Lagrange equations 
(Bl) one obtains the equations of motion (10) written 
in terms of the collective variables. 

Let us now consider energy of the system £. Substi¬ 
tuting the Ansatz (7) into the energy expression (4) and 
performing the integration one obtains (up to an additive 


constant) 



£ e 



^ “ A 


pFi cos $-sin $ 

2 

OO 


^ ^ cosh ^ 


r(s)ds. 


— OO 


Static form of Eqs. (3) reads 6£./69 = 0, SS./S(j) = 0. 
Taking into account form of energy (4) and exchange en¬ 
ergy density (5) one obtains the following set of equations 


6" — sin 9 cos 9 sin^ 4> + w (Ala) 

— 2())'k sin^ 9 + k' cos (j) = 0, 

+29' [cot 9(j)' + K sin (j)] + sin </> cos (j) (Alb) 
—k' cot 9 sin (j) = 0. 

Equation (Alb) has solutions (j) = cj)o = 0, tt. Substitu¬ 
tion of this solution into (Ala) results in equation (6), 
which determines the function 9{s). 


Appendix B: Details of the g — <!> model usage 

The equations of motion (3) are the Euler-Lagrange 
equations 


<5/; _ d _ (5T 

dt sii si ’ 




(Bl) 


In the limit case kA 0 the expression (B5) is re¬ 
duced to (8). For the more general case 0 < Ak < 1 
one has / < 1 and consequently Fi > F 2 . In this case 
the value of phase $0 which minimizes the energy (B5) is 
determined as cos $0 = —p and the corresponding value 
of the equilibrium domain wall position qq is determined 
by the equation (15). The case Ak > 1 (domain wall 
width is larger than the curvature radius) is not consid¬ 
ered here, because the application of the Ansatz (7) is 
questionable in this case and a precise solution of equa¬ 
tion (6) is required. 

The equations of motion linearized in vicinity of the 
equilibrium q = qo, 4) = $0 read 


(l + a^) 


q 

!> 



aA 

A 

Ao 


Ao® 

—o3 


q 

!> ’ 


(B6) 


where q = q—qo and 4) = 4)—4>o are small deviations from 
the equilibrium and quantities A and ® are defined in 
(14). For the case of small a the solution of (B6) results in 
harmonic decaying oscillations q oc sin(Ut-|-(5o)e“''*, 4> oc 
cos{nt + 6o)e~'^* with frequency fl and effective friction 
r] presented in (14). 
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For the case of parabolic wire (16) one obtains A = 


yi(AoKo) and ® = ®(AoKo), where 




%{x) = 


1 — X 


sech 


/(«) 


(l+^2)i 1+^2 


sech^f^ 


(B7) 


-de 


These expressions are used to calculate the eigenfrequen- 
cies and friction parameters in (14), see solid lines in 
Fig. 3. 
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